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ABSTRACT 

Context. The spectral features of HD 94509 are highly unusual, adding an extreme to the zoo of Be and shell stars. The shell 
dominates the spectrum, showing lines typical for spectral types mid-A to early-F, while the presence of a late/mid B-type central star 
is indicated by photospheric hydrogen line wings and helium lines. Numerous metallic absorption lines have hroad wings hut taper to 
narrow cores. They cannot be fit by Voigt profiles. 

Aims. To describe and illustrate unusual spectral features of this star, and make rough calculations to estimate physical conditions and 
abundances in the shell. Furthermore, the central star is characterized. 

Methods. We assume mean conditions for the shell. An electron density estimate is made from the Inglis-Teller formula. Excitation 
temperatures and column densities for Fe i and Fe ii are derived from curves of growth. The neutral H column density is estimated 
from high Paschen members. The column densities are compared with calculations made with the photoionization code Cloudy. 
Atmospheric parameters of the central star are constrained employing non-LTE spectrum synthesis. 

Results. Overall chemical abundances are close to solar. Column densities of the dominant ions of several elements, as well as 
excitation temperatures and the mean electron density are well accounted for by a simple model. Several features, including tbe 
degree of ionization, are less well described. 

Conclusions. HD 94509 is a Be star with a stable shell, close to the terminal-age main sequence. The dynamical state of the shell and 
the unusually shaped, but symmetric line profiles, require a separate study. 
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1. Introduction 

3 

The unusual nature of the spectrum of HD 94509 was noted 
during a search of the ESO archive for possib le Herbig Ae/Be 
stars. S imbad lists HD 94509 as a Be star, while lHouk & Cowled 
lll975l) call it Bp SHELL, with the remark that it has strong typ¬ 
ical shell lines, with Hj3 somewhat filled in. The star appears 
among a numbe r of surveys for young or unusual objects, e.g. 
The et alj (Il994t) and iReed & Beatt'^ (Il995h . that give little spec¬ 
troscopic description beyond a basic classification. 

_ H D 94509 is a shell st ar by the modern criteria dHanuschikl 

Il996t iRivinius et ani2006l) . which required the central absorp¬ 
tion of the emission lines to be below the level of the underlying 
stellar continuum. The typical picture of a shell star discussed in 
these papers was with a flattened shell close to the line of sight 
or nearly 90° to the plane of the sky. 

IWade et ^ (l2007t) included HD 94509 among their survey 
of possible Herbig AeBe stars for magnetic fields, but their pa¬ 
per does not comment on ot her characteristics of th e spectrum. 
It appears in the catalogue of ide Winter et alJ (l200ll) as a Herbig 
Be star. However, the raw infrared colors (J - H - -I-0.054, and 
H - K- -0.061) place HD 94509 within the group of classical 
Be stars, and well displaced from the region of Herbig Ae/Be 


* Based in part on ESO programs 082.C-0831, 083.A-9013, and 
084.C-0952 


stars delineated bv iHernandez et akl (l2005h . Reddening correc¬ 
tions would only make the displacement larger. But dilution in¬ 
dicators (see below) show that HD 94509 is a shell star rather 
than a classical Be star. 

We could And no detail ed study of the spectrum, though 
ICorporon & Lagrange! (1 19991 henceforth, CL) discuss HD 94509 
in a note on individual sources in their survey of binaries among 
Herbig Ae/Be stars. They noted He 1/16678 as well as neu¬ 
tral lines of Ca i and Le i, and called the spectrum composite. 
We confirm the underlying He i line, but see no evidence of a 
metallic-line spectrum from a secondary star. CL did not discuss 
the metallic-line profiles we And so remarkable. 


2. A puzzling spectrum 

2.1. Observations 

The study is based primarily on UVES archive spectra of HD 
94509 obtained with the dichroic beam splitter DIC-2 in blue 
and red arms on 16 March 2009. The high-resolution spec- 
trograph UVES (U ltraviolet and Visual Echelle Spectrograph, 
iDelAer et al.l 120061) is installed at the VLT-8m (ESO/Paranal, 
Chile). Our wavelength coverage is 3750—4985, and 5682- 
9460 A. The resolving power is/? = A/AA 40000, due to 
the use of a wide slit. The signal-to-noise (S/N) varies but was 
typically about 330 in the blue and 250 in the visual. 
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Fig. 1. Comparison of the spectrum of HD 94509 (lower curve) with a FEROS spectrum of o Aqr (B7-B8 IV-V, upper curve), a prototype Be shell 
star, for a selection of strategic lines. The most prominent features are identified. Both spectra were rv-corrected to the laboratory restframe. 


We also used an ESO archive FEROS spectrum (R - 48 000, 
S /N ^ 120) obtained on 7 June 2009, covering the wavelength 
range from 3666 to 87 00 A. FEROS (Fibe r-fed Extended Range 
Optical Spectrograph. iKaufer et al.lll999l) is a high-resolution, 
Echelle spectrograph operated at ESO in La Silla, Chile and is 
installed at the MPIA 2.2m telescope. Both UVES and FEROS 
spectra were examined for spectral variations. The FEROS spec¬ 
trum was also useful for features (e.g. Ban 44934) that were 
in UVES order gaps. Additionaly, an X-Shooter spectrum of 
HD 94509 from 5 February 2010 was retrieved from the ESO 
archive. X-Shooter (I Vernet et al.ll20 111) is installed at VLT-8m 
(ESO/Paranal, Chile) and is the first of the second generation in¬ 
struments at VET. The resolution and noise level of this spectrum 
was considerably lower than those of the UVES and FEROS 
spectra. We only used the X-Shooter data to look for spectral 
variations. 

Figure [1] shows a comparison of the spectrum of HD 94509 
with that of o Aqr (HD 209409, B7-B8 V-IV), a prototype Be 
shell star, for a selection of strategic lines. While some similar¬ 
ities in the Balmer and He i lines are seen, the richness of the 
metallic-line spectrum and the peculiar line profile shapes stand 
out, and - to our knowledge - are unreported for any other star 
so far. Note that o Aqr shows sharp line cores originating from 
the shell in a few of the strongest Fe ii lines, e.g. 44924. 

CL discuss three absorption lines in the neighborhood of the 
broad Hei line 46678, which they tentatively identify as Fei 
46677.989, Cai 46690.628, and Cai 46717.681, with velocity 
shifts of -45, -35, and -7 kms^'. Their Fig. 20 shows wave¬ 
length shifts at three epochs, best seen in the strongest feature. 



Fig. 2. UVES (black) and FEROS (gray/red online) spectra in the re¬ 
gion of CD’s Eig. 20. Three absorption features are identified as Eeii 
(46677.30), Tin (46680.13), and Tin (46717.79). 

The upper part of CL’s Fig. 20 closely resembles our spectra, 
shown in detail in Fig. |2] The broad He i absorption is appar¬ 
ent. We can plausibly identify the sharper absorptions as lines 
of Fe II and Ti ii with the radial velocity of the other metallic 
lines of the shell. Interestingly, 46677 of Fe ii does not appear 
in the current National Institute of Standards and Technology 
(NISTB list, though it is in the Muli^let Tables, and in the Vi¬ 
enna Atomic Line Database (VALDfl G. Nave of NIST (private 

^ http://www.nist.gov/pml/data/asd.cfm 
^ http://vald.astro.univie.ac.at/~vald/php/vald.php 
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Fig. 3. The strong Feii line .14233, illustrating the unusual metallic-line 
profiles. The gray (red online) synthesis is based on a = 10000 K, 
log g = 2.5 non-LTE spectrum synthesis model with 10 x solar abun¬ 
dances, but with the Stark broadening increased by a factor of 6 x 10^! 
The dark gray (blue online) rotational profile used a 10 x solar iron 
abundance, with V ■ sin(i) =100 km s“' and standard Stark broadening. 



4220 4225 4230 4235 4240 4245 4250 

Wavelength [A] 


Fig. 4. To illustrate the richness of the metallic spectrum of HR 94509, 
we compare the region of the strong Feii line 44233, in HD 94509 
(black) and that of the F3 la supergiant HD 74180 (gray, red online). 
The somewhat rectangular shape of the profiles with rounded bottoms, 
of the latter star may be rather closely matched by a model with appro¬ 
priate microturbulent and rotational broadening. 


communication) suggests it is masked in her spectra by an argon 
line. Our material does not indicate the presence of another star, 
or another shell with a relative radial velocity shift. We can, of 
course, not exclude a stage of variability at an epoch not covered 
by our spectra. 

2.2. The metallic lines 

Figure[3is typical of the stronger metallic-line profiles. A rough 
fit to the observed line is possible with a 10 x solar abundance 
and a standard model atmosphere see Sect. [3 for details on the 
modelling. However, the wings require an impossibly large en¬ 
hancement of the damping, by 6 x 10^. The figure also shows a 
calculated profile based on a V • sin(/) of 100km s 'and standard 
Stark broadening. 

The stronger metallic-line profiles of HD 94509 gracefully 
taper to cores as sharp as the instrumental resolution allows. 
The profiles merge with the continuum in what resemble Lorentz 
wings. The profiles are very symmetrical, with little or no indi¬ 
cation of a wind, or velocity asymmetry within the zone in which 
they are formed. 

A broader view of the /14233-region is shown in Fig. |4l The 
richness of the HD 94509 spectrum is comparable to that of 
the F3 la supergiant HD 74180, and not common among shell 
stars, cf. Fig. [1] Note we do not imply that HD 94509 is a com¬ 
posite containing such a star. Dilution indicators, discussed be¬ 
low (Section [4. II) . show that the metallic-line spectru m of HD 
94509 is not photospheric. But a f ew stabl e shell stars (iGulliveii 
Il98lh may have similar spectra. iMerrilll (Il952h shows a mi¬ 
crophotometer tracing for HD 193182, whose spectrum is surely 
rich in metal lines. We have been unable to find a digital, high- 
resolution spectrum that would allow us to make a close com¬ 
parison of the line shapes. 

A few of the strong metallic lines show broad emission 
shoulders as well as relatively narrow, deep absorption. An ex- 
ample is the F e ii emission, /16516: a®S 2 . 5 -z®D 3 ^ in Multiple! 40 
(lMoorelll945h . Multiple! 42 of Feii (TT4924, 5018, and 5169) 
is commonly in emission in many peculiar spectra (Be, Herbig 

Note that the effective temperature T^g and surface gravity logg of 
the model are chosen here for illustrative purposes only - the line is not 
formed in the atmosphere of the central star, but in the shell. 
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Fig. 5. Base, or full width of 31 Feii absorption lines. An average half 
width for these lines corresponds to a velocity of some lOOkms^'. No 
correlation with wavelength is seen. See text for uncertainties in mea¬ 
surements 

Ae/Be, T Tauri). It is is also a sextet to sextet transition involving 
the same lower (a®S) term. Interestingly, Multiple! 42 is at most 
only weakly in emission in HD 94509. 

In this paper, we will not attempt to account for the overall 
line shapes, especially the graceful tapering toward the cores of 
the stronger metallic lines. Figure [2 shows that velocities of the 
order of 100 km s 'are required to account for the line wings. If 
the Doppler effect were responsible for the width of these lines, 
a correlation of the width with wavelength would be expected. 
However, such a correlation is not obvious in Fig. |5] 

Measurements of the base widths of absorption lines are 
most uncertain, because a qualitative judgement must be made 
of the position where the line merges with the continuum. The 
difficulty is compounded by the presence of blending. Remea¬ 
surements of a sample of 10 randomly chosen Fe ii lines shows 
an average difference of 12%, with an extreme value of 23%. In 
spite of the uncertainty of these base widths, we do find a signif¬ 
icant correlation of the base width with the line depth, as shown 
inFig.|6] 

We offer no explanation of the unusual line profiles or the 
correlation of Fig. |6] Widths due to turbulent macroscopic mo¬ 
tions would require Mach numbers of 8 or more, and seem un- 
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Fig. 6. The base, or full width of Fe ii absorption lines shows a definite 
correlation with the line depths. For 31 points, the Pearson correlation 
coefficient is 0.674; the significance, or probability of no correlation is 
3 X 10-^ 



Fig. 7. The Paschen line, P12, shows the sharp core of the shell and 
the broad-line morphology of the central source. Overlaid are rotational 
profiles for 0 (dark gray, or blue) and 260 km s“* (light gray, or red). Two 
NI lines are seen to the violet of P12. See the text for details. 


likely. Rotation of the central star, or Keplerian revolution of 
a disk could readily produce widths comparable to those mea¬ 
sured. It is puzzling that we cannot establish a clear correlation 
of the line widths with wavelength. 

Possibly calc u lations such as those performed by 
iHummel & Dachs] (Il992h . but for absorption lines, could 
account for the line shapes. We must leave this for future work. 

2.3. The hydrogen lines 

The B aimer lines are highly symmetrical. Both Ha and and H/3 
have double emission peaks of equal intensity, see Fig. [T] The 
central absorption cores yield velocities that closely agree with 
the metallic-line absorption. The central reversals of both Ha 
and Hy6 are very deep, while the cores of Hy, through He are vir¬ 
tually black. Our spectra extend from Ha to Ht or HIO, and none 
of the core wavelengths deviate from the laboratory positions 
by more than 0.04 A. For such broad profiles, this may easily 
be attributed to measuring error. The mean deviation for the 9 
lower Balmer cores from that expected from the metallic cores 
is 0.0033 A. The symmetry of both the hydrogen and metallic 
absorption cores is indicative o f “negligible line-of- sight veloci¬ 
ties” in the absorbing shell (cf. fevinius et al.1l2006l) . 
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Fig. 8. Like Fig. |7] but for Hy. The black core of the line is indicative 
for the shell covering the entire stellar disk. 

Table 1. Apparent rotational velocities (half widths) in kms^'. Emis¬ 
sion lines marked with an asterisk. An average is given for the Can 
infrared triplet. 


Spectrum 

Wavelength (A) 

V ■ sin(i) 

He I 

5876 

255 

He I 

6678 

299 

Sin 

6347, 6371 

245 

Caii(IRT)* 

8500-8660 

225 

Ave. 31 Fell 

3783-7516 

99±33sd 

FeiT 

6516 

214 

Ki 

7699 

114 


Notes. * The width measurement is of emission shoulders. 


The non-LTE calculations of PI2 shown in Fig. [T] were made 
with a 14000K, logg = 3.5 model,assuming an additional con¬ 
tinuum equal to a factor of 0.2 times the stellar continuum. The¬ 
oretical plots are for V ■ sin(/) = 0 and 260 km s“’. As the Stark- 
broadened wings of P12 are of stellar origin (only the sharp line 
core originates in the shell), this may be viewed in support of the 
atmospheric parameters chosen for the central star (Sect. lO l. 
The same applies to other hydrogen lines, also the Balmer lines 
from Hy (see Fig.O to the higher series members. A continuum 
contribution of 0.12 times the central star’s photospheric contin¬ 
uum was adopted for Hy, meeting the expectation of a decreas¬ 
ing ratio towards the blue. Note, however, that the true amount 
of dilution of the underlying spectrum is uncertain. 

2.4. The helium lines 

The He i lines are broad and shallow, unlike the metal lines. They 
can be matched by rotationally-broadend line profiles, except 
for cases with metal line blends (see e.g. Fig. |2]l. This is indi¬ 
cation for a stellar atmospheric origin of these lines. They can 
be fit using atmospheric parameters as adopted for the central 
star (Sect. lOT i. and a solar helium abundance. Figure |9] shows 
the example of the line Hei 45875, assuming an additional con¬ 
tinuum from the shell equal to a factor of 0.15 times the stellar 
continuum, intermediate between the two values adopted for the 
hydrogen Balmer and Paschen line above. A rotational velocity 
of ~260kms ' is indicated, in consistency with the typical pic¬ 
ture of shell star, where the metallic absorption is superimposed 
on the spectrum of a rapidly rotating central star. 
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Fig. 9. The He i /15875 line, which originates in the atmosphere of the 
central star, with no contribution from the shell. A rotational profile 
for 260kms“* (light gray, or red) fits the observation. See the text for 
details. 

2.5. Additional remarks 

Estimates of the putative V ■ sin(/) may be made from the (half) 
base widths of both additional emission and absorption lines. 
Uncertainties in these values depend on the line width measure¬ 
ments, whose accuracies are discussed in Section lZ^ 

While the metallic spectrum is predominantly absorption, a 
few lines show “emission shoulders”. Table [1] summarizes the 
measurements, yielding the following picture. The lines origi¬ 
nating from the highest excitation potentials (i.e. the Hei lines) 
indicate the highest apparent rotational velocities, of about 250- 
300 km s '. This should correspond to the rotational velocity of 
the central star, which is in line with a typical V ■ sin(/) that 
range s from 200 to 450 km s ' for the central component of shell 
stars dRivinius et al.l20()5) . Slightly lower velocities are obtained 
from the emission shoulders of spectral lines, which arise from 
the faster-spinning inner parts of the disk. Then, even lower ro¬ 
tational velocities are found for the metallic absorption lines, 
which involve low-excitation levels and indicate the formation 
in the slower-rotating outer parts of the disk. 

Finally, we note that several broad and shallow features 
are seen in the spectrum. These coincide in wavelength with 
several of the typically most pronounced diffuse interstellar 
bands (DIBs). In particular /1/15780.45, 6283.86 and 6613.62 are 
present, but shifted by ~ H-17 km s“' relative to the intrinsic spec¬ 
trum of HD 94509. This is is consistent with the findings from 
the sodium D and calcium H and K lines (see Sect. 0. 

3. Characterization of the centrai star 

While the overall spectrum is dominated by shell absorption and 
emission, a small number of features are of photospheric ori¬ 
gin, as indicated in the last section: the wings of several hydro¬ 
gen lines and Hei lines. Photospheric metal lines unaffected by 
disk absorption cannot be identified unambiguously. These lines 
do not provide sufficient information for a fully-fledged quan¬ 
titative analysis, but allow an estimate of the atmospheric pa¬ 
rameters. A Tgg- 14 000 K and logg = 3.5, solar helium abun¬ 
dance and V ■ sin(/)»260kms ' seems a proper choice, but 
cannot be considered precise. The presence of Hei lines set a 
lower boundary of the temperature to ~ 11 000 K. However, this 
would require an enhanced helium abundance, which is unlikely 
in view of evolution models for rotating stars that do not pre¬ 
dict a noticeable increase of the surface helium values even at 



Fig. 10. Dilution in the shell of HD 94509. The Sin lines 4/14128 and 
4130 (arrows) are stronger in the A2 lae star Deneb (grey, red on¬ 
line) than in HD 94509 (black) because of dilution. The absorption at 
/t4128.12 in HD 94509 is largely due to /14128.13 of Mnii. 

critical rotation in the relevant mass range (iBrott et al.ll2011l : 
iGeorgy et al.l2013l) . The wings of the hydrogen lines provide ad¬ 
ditional information on the effective temperature, but are subject 
to temperature-gravity degeneracy. Calculations with the ioniza¬ 
tion code CLOUDY (see Section |6]l show that the Ha emission 
strength would be too large for a higher temperature than the one 
adopted. 

The adopted gravity value implies that the star is in an ad¬ 
vanced stage of the main-seq uence evolution, con sistent with 
theoretical considerations (e.g. iGranada et al]l2013h . Note that 
the light from the central star - as a consequence of the fast 
rotation and the shadowing by the disk seen equator-on (i.e. 
i ~ 90°) - is dominated by the hotter and slower-rotating po¬ 
lar caps. Consequently, our above-mentioned parameters are 
likely upper values for the non-uniform Tgff- and logg- distri¬ 
bution over the rotationally-flattened ste llar surface, while the 
rotational velocity is a lower limit (e.g. iTownsend et al.l 120041 : 
iFremat et al.l2005h . The metallicity of the star seems close to so¬ 
lar, o r compatible with the pre sent-day cosmic abundance stan¬ 
dard (iNieva & Przvbill^l2012l) . as indicated by the analysis of 
the shell (see Sect.|6]). 

We employed a hybrid non-LTE approach for our stellar pho¬ 
tospheric line-formation calculations, which has been success¬ 
fully use d before in a similar temperature regime to the o ne stud¬ 
ied here dPrzvbilla et al.ll2006l : iNieva & Przvbillall200'^ . Plane- 
parallel, hydrostatic and homogeneous L TE model atmo spheres 
were computed with the code Atlas9 (iKuruczl 1199Ih . Then, 
non-LTE line-formation co mputations were performed with the 
codes Detail and Surface (lGiddingslll98ll : Butler & Gidding j 
119851 both updated and improved). Detail calculates atomic- 
level populations by solving the coupled radiative transfer and 
statistical equilibrium equations, and Surface computes the for¬ 
mal solution using rea listic line-broadening func tions. We em- 
pl oyed model atom s of iPrzvbi lla & jut^ ([gOO^ for hydrogen, 
of lErzybilial (l2005h for He, and of Becked dl998h for Fe ii. 

We note that the model calculations are intended to illustrate 
the plausibility of our parameter choices, but are well aware that 
this approach is suited to address the complex phenomenology 
of a Be (shell) star only coarsly. We also note that similar ap¬ 
proaches are nevertheless employed for quantitative Be star anal- 
yses (e.g. lLevenhagen & Leisteil2004t iDunstall et al.ll201 ih . 
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Table 2. Electron densities log(A^j) 


constant 

n = 41 

n = 42 

23.26 

11.16 

11.08 

23.49 

11.40 

11.32 


4. Physical conditions in the shell 

4.1. Dilution indicators 

A general characteristic of shell stars is the relative enhancement 
of second spectra (from ions) relative to first (from neutrals), 
when compared with the spectra of dwarfs or normal giants. This 
may be attributed to the low electron pressure in the shell, which 
will be discussed in more detail below (cf. Sect. I4.2b . A sec¬ 
ond charact eristic, discus sed in numerous papers and reviews by 
Struve (e.g. lStru^ll95ll) . is the weakness of certain lines which 
he called “dilution indicators”. Specifically, Struve called atten¬ 
tion to lines of Siii /1/13853-3862, 4128-4130, and Mgii 4481. 
He noted that the lower levels of the dilution indicators were 
connected to the ground level by strong permitted transitions. 
Most absorption lines of Fe ii and Ti ii arise from levels of the 
same parity as the ground term, and are thus metastable. 

We illustrate this in Fig. [TO] where lines in HD 94509 and 
the supergia nt Deneb are shown. The la tter spectrum is the one 
described bv ISchiller & Przvbillal(l2008l) . Generally, the metallic 
lines in HD 94509 are significantly stronger than in Deneb, but 
this is not the case for Struve’s dilution indicators. 

4.2. Electron density 

A depth-dependent model of the shell of HD 94509 is beyond 
the scope of the present work. We attempt to determine mean 
conditions for the shell as a whole, making use of techniques 
common for normal stellar atmospheres prior to the use of nu¬ 
merical modeling. 

Our high-resolution spectra do not include the Balmer limit, 
but the Paschen limit is clearly seen. The line P41 was measured, 
but P42 was not, and an inspection of the spectrum reveals no ob¬ 
vious feature between P41 and three moderate N i lines. A well- 
established interpretation of the termination of a one-electron 
line series is that overlappin g Stark wings merge s o individual 
features are no longer seen (linglis & Telle^ Il939l he nceforth, 
IT). A simple derivation of the formula may be found in lCowle^ 
(Il970l) . where it is assumed that protons and electrons contribute 
equally to the line merging, unlike IT. If the last Paschen line is 
Pn, then 

log(A,) = log(Ap) - 23.26 - 7.5 log(n), (1) 

where Ng an d Nn are the el ectron and proton number densities, 
respectively. I Varslinil (1 1 978h uses the constant 23.49 rather than 
23.26. Let us put n - 41 and 42 into Eq.[T] Results are shown in 
Table |2] We conclude that a mean logarithmic electron density 
of the shell is 11.2+0.5. 

4.3. The shell temperature 
4.3.1. Fell 

We know from the dilution indicators that there can be no sin¬ 
gle temperature that will describe the excitation and ioniza¬ 
tion within the shell. Nevertheless, the ionization and excita¬ 
tion may be approximately estimated from measurements of the 



Fig. 11. Comparison of ILi’s (mA) measured using segmented profile 
fits and those using fits to a Voigt profile, the average absolute value of 
the differences in the plotted equivalent widths is 14 mA 


Table 3. Excitation temperatures. See text for uncertainty estimates. 


Spectrum 

No. lines 

^(eV) range 

Tex (K) 

Uncertainty (K) 

Tin 

26 

0.62-3.12 

9600 

600 

Fei 

29 

0.00 - 3.63 

7750 

650 

Fell 

46 

2.58 - 7.27 

9164 

170 


line strengths. These estimates will be subject to clarification by 
models of the shell. We have worked primarily with the second 
spectra of metallic lines; Sc ii, Ti ii. Cm, Fe i, Fe ii, and Ni ii. The 
Fe II spectrum is by far best suited for our purposes, both because 
of its richness and the detail of its treatmen t by the versatile pro¬ 
gram Cloudy (cl3.03. lFerland et al.ll2013l cf. Sect. 16.11 1. 

The unusual shape of the Fe ii (and other metallic line) pro¬ 
files means that measurements of equivalent widths Wa require 
special care. Our general approach in abundance work is to fit 
metal lines with Voigt profiles. However, the a - yUlAAo) is 
usually small, typically 0.01 to 0.1, so the profiles are not far 
from Gaussian. Such a low a-value is completely unsuitable for 
HD 94509. For the lines with profiles that are not too deep, we 
could make tolerable fits using a - 5, virtually Lorentzian pro¬ 
files. Core and wings of the stronger (deeper) lines could not 
be fit by any Voigt profile. All of the profiles were also fit with 
small line segments, so for many lines we could compare the 
results (Fig.iriTi. We adopted simple means for all pairs of mea¬ 
surements. 

One may estimate an overall excitation temperature by plot¬ 
ting log{WA/A) vs. logigfA) - (5040/r) ■ where is the 
lower excitation potential of the absorption line. One adopts 
the temperature T as the value that gives the minimum sepa¬ 
ration of lines with different x- We illustrate this determ i nation 
using reliable Fe ii g/-values, from iMelendez & Barbil^ (l2009l 
cf. also Adelman 2014). Figure [12] shows the adopted fit, using 
6 - 5040/7’ = 0.55, or Texe = 9164 + 170 K. The uncertainty 
is derived from trial plots with 6 = 0.55 and 0.56, which were 
judged marginally acceptable. 

4.3.2. Fei and Ti ii 

Other atomic species with enough measurable (mostly un¬ 
blended) lines to allow a credible curve of growth are Fe i and 
Ti II. Results are gathered in Table |2 uncertainty estimates, ex¬ 
plained in 14. 3. II are significantly larger for Ti ii and Fe i, which 
have fewer lines than Fe ii. 
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Log(gfA)-0x 

Fig. 12. Excitation temperature from Feii lines in the shell. Colors may 
be seen online. The (red) stars indicate lines with the highest lower exci¬ 
tation energies (6.2 to 7.3 eV). The lowest-excitation lines (2.6-2.9 eV) 
are indicated by black squares. The black line indicates a 45° slope. 
This plot was made with 9 = 0.55 or T = 9164 K with an uncertainty 
of ~170K (see text). 

4.4. Column densities 

The column density of an absorber follows from the relation for 
weak lines with equivalent widths W'. 

2 

IT' - ^A^Nnfnm- (2) 

mc^ 

Here, N„ is the number (per cm“^, column density) in the lower 
energy level, in the absorbing shell. We have assumed the metal 
line absorption arises entirely in the shell. In this context we note 
that the central star has a spectral type in the late/mid-B range 
and it is unlikely to find chemical peculiarities, as the fast rota¬ 
tion inhibits the CP phenomenon because of rotationally-driven 
hydrodynamical mixing. Therefore, one may expect the only fea¬ 
tures of iron-group species present in the stellar spectrum to be 
weak Fe ii lines. At most, these stellar lines contribute a few per¬ 
cent to the measured equivalent widths, i.e. they are indeed neg¬ 
ligible. In fact the presence of strong unusually-shaped features 
of e.g. Sen, Tin, Cm, Fei and Niii makes the spectrum appear 
so unusual. 

4.4.1. Fell 

One may estimate from Fig. [12] that lines as strong as 250 mA 
(\ogiWxlA + 6 = 1.7, for 5000 A) are still on the linear part of the 
curve of growth. Let us make the provisional assumption that 
iron is predominantly first ionized, neglecting Fe i and Fe iii. 

Using the Boltzmann formula, with u{T) the partition func¬ 
tion of Fe^, and using logarithms, this may be written to give the 
column density of all Fe^ for each line. With 0 = 5040/r, and 
the lower excitation energy ;t'„, we have 

log[A(Fe+)] = log(lT^)//l)-Hl2.O53-log(g/T)-Hlog(M(7’))H-0^„. 

(3) 

We work out the column densities for Fe^ using Eq.j^for the 
Fe II lines shown in Tab.|4| The average value (of the logarithms) 
of the last column is 18.77+0.19. 

4.4.2. Neutral iron 

While the Fe i lines are less numerous than those of Fe ii, a num¬ 
ber of them may be measured, and analyzed in a similar manner. 


Table 4. Calculation of column densities for Fe+, Dimensions are given 
for the arguments of the logarithmic quantities. 


A 

A 

T 

eV 

log(g/) 

no units 

logiVf,] 

mA 

log[A(Fe+)] 

cm“^ 

4147.27 

4.616 

-3.508 

0.88 

18.54 

4413.60 

2.676 

-3.79 

1.76 

18.58 

4601.40 

2.891 

-4.48 

1.34 

18.93 

4648.94 

2.583 

-4.58 

1.54 

19.05 

4720.15 

3.197 

-4.48 

1.22 

18.96 

4833.20 

2.657 

-4.64 

1.31 

18.89 

4840.00 

2.676 

-4.75 

1.25 

18.95 

4893.82 

2.828 

-4.21 

1.67 

18.90 

5813.67 

5.571 

-2.51 

1.65 

18.54 

6084.11 

3.199 

-3.79 

2.10 

18.93 

6103.51 

6.217 

-2.171 

1.93 

18.80 

6317.40 

6.223 

-2.16 

1.84 

18.67 

6627.22 

7.274 

-1.609 

2.09 

18.91 

7301.57 

3.892 

-3.63 

1.85 

18.74 

7515.81 

3.903 

-3.39 

1.99 

18.62 

7841.40 

3.903 

-3.721 

1.42 

18.35 

average 




18.77+0.19 (sd) 


Table 5. Calculation of column densities for neutral Fe or Fe”. The for¬ 
mat is the same as for Table|4] 


A 

A 

T 

eV 

log(g/) 
no units 

logiVf,] 

mA 

log[A(Fe°)] 

cm“^ 

3873.76 

2.43 

-0.88 

1.72 

15.52 

3983.96 

2.73 

-1.02 

1.45 

15.54 

4118.55 

3.63 

0.22 

1.88 

15.20 

4134.68 

2.83 

-0.68 

1.53 

15.30 

4175.64 

2.84 

-0.83 

1.32 

15.24 

4181.76 

2.83 

0.37 

1.75 

14.46 

4309.37 

2.95 

-1.19 

1.83 

16.15 

4442.34 

2.20 

-1.26 

1.18 

15.11 

4459.12 

2.18 

-1.28 

1.48 

15.42 

4466.55 

2.83 

-0.60 

1.69 

15.32 

4476.02 

2.84 

-0.82 

1.52 

15.37 

4494.56 

2.20 

-1.14 

1.62 

15.43 

8688.63 

2.18 

-1.21 

1.96 

15.26 

average 




15.33+0.35 (sd) 


Oscillator strengths for Fe i are from iFuhr & Wies^ (l2006l) . Ta¬ 
ble |5] presents the results. The average of the logarithms for the 
last column is 15.41+0.35. 


4.4.3. Sc II, Ti II, Cr ii, and Ni ii 

Column densities have been estimated for other species, using 
the Fe ii curve of growth to choose “weak”’ lines or estimate sat¬ 
uration corrections. Results are discussed in Sect.|6|below. These 
are uncertain by at least a factor of three because of the difficulty 
in estimating the contributions from the wings, especially in the 
case where wings from different lines overlap. Non-LTE effects 
are also expected, as can be seen from the more detailed results 
for Ee II available from Cloudy (cf. Sect.|6]l. 

A resume of observed and theoretical column densities may 
be found in Sect. |6] Relevant data are presented in Table lA.il in 
Appendix|A| 
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Table 6. Calculation of column densities (cm for « = 3, from 
Paschen lines, P«, where n is in the first column. 


n 

WAmA) 

/ 

loglAj) 

30 

368. 

2.034E-04 

15.474 

31 

354. 

1.840E-04 

15.501 

32 

259. 

1.670E-04 

15.408 

33 

210. 

1.521E-04 

15.358 

34 

186. 

1.388E-04 

15.345 

35 


1.271E-04 


36 

160. 

1.167E-04 

15.355 

37 

138. 

1.075E-04 

15.327 

38 

116. 

9.901E-05 

15.288 

39 

90. 

9.150E-05 

15.212 

40 

89. 

8.474E-05 

15.241 

41 

73. 

7.863E-05 

15.188 

42 

61. 

7.309E-05 

15.138 

average 


15.32±0.11 (sd) 


4.4.4. Hydrogen 

The highest Paschen lines (those that have no contribution from 
the central star) may be used to make a crude estimate of the col¬ 
umn density N^, of hydrogen in the n -3 level. Measurements of 
equivalent widths were made using the segmented fits described 
above. With Eq.|2] we obtain values of Nj, that would apply as 
if the lines wer e unsaturated. In the ideal case, discussed in the 
classical text of lUnsoldl (Il955h . we would expect values of to 
increase, as the lines become optically thin, and then to decrease 
as the lines blend together, so their entire width is not measured. 
Results are shown in Table |6] The line P35 was blended and 
hence not used. Absorption oscillator strengths were calculated 
using exact expressions for the wave functions and relevant inte¬ 
grals. 

The last column of Table |6] decreases as n increases, as ex¬ 
pected. We do not see a decrease due to saturation for the smaller 
n. However, the variation from P30 to the last visible line, P42, 
is not large, so the uncertainty is comparable to other approxi¬ 
mations made herein. We adopt log(A 3 ) = 15.3. 

5. Qualitative chemistry of the circumstellar 
envelope or shell 

Of the lighter elements, hydrogen is well represented by the fea¬ 
tures described previously. We see broad absorption from He i 
/15876 and 6678, and most of the bluer optical lines, but this is 
from the underlying stellar spectrum. There is no sign of emis¬ 
sion in the He i features. There is also not a hint of absorption or 
emission at the position of the Li i resonance lines near 46708 - 
which is expected because any lithium initially present would 
be transported to deeper layers of the stellar envelope due to 
rotationally-induced mixing, and destroyed by nuclear reactions. 
This applies in a similar manner to beryllium and boron (relevant 
lines are located outside the available spectra). Absorption near 
/1/19062 and 9112 could be due to C i in Multiplet 3. No C ii lines 
were noted. A number of absorptions due to N i were identified 
in the near infrared. Several of these appear in Figs.[T]and|7l The 
sharp cores of the N i lines are indicative of absorption from the 
shell. Neutral oxygen lines from the shell are the strong infrared 
triplet 447772, 7774, and 7775 and the line 48446 (see Fig.[B. 

The D lines of sodium are strong. Each line is double, with 
the deeper component having the stellar radial velocity, with 
the strong, second component shifted to the red by 17.3 km s '. 


Table 7. Characteristics of the Cloudy model which fit the observations 
as summarized in Table[8l Keplerian velocities, Vk are in kms^', for a 
4.2 Mg star. 


Central star 

T,g = 14000K, logg = 3.5 

Stellar luminosity 

2.4 X 10^® erg s^' 

Stellar radius 

4.2 Rg 

Inner cloud radius 

14.4Rg; Vk = 236 km s-' 

Outer cloud radius 

45.5 Rg; Vk = 133 km s-' 

Cylinder semi-height 

5.0 Rg 

Total hydrogen density 

1.58 • 10" cm“^ (constant) 

Average electron density 

9.8 • 10*° cm“^ (variable) 

Molecules and grains 

off 


likely being of interstellar origin (an analogous situation is seen 
in the calcium H and K lines). Several lines of Mgi are identi¬ 
fied including the b-lines, which have the broad wings and sharp 
cores of the shell. The usually strong Mgii 44481 is severely 
weakened as noted earlier. The A1 1 resonance doublet is present, 
but weak. No Si i lines were identified, but Si ii lines are present 
in addition to the dilution indicators mentioned in Sect. 14. II The 
violet component of the Ki resonance pair, 47665 is obscured 
by water-vapor absorption, but 47699 is unblended and strong 
(Wi ~500niA), see Fig. [T] The 3d elements Ca-Ni are all well 
represented, but heavier elements are weak or absent altogether, 
similar to the situation in a supergiant spectrum at comparable 
temperature. 

An analysis by wavelength coincidence statistics (WCS) re¬ 
veals meager evidence for the presence of neutron addition ele¬ 
ments. Even the Sr ii resonance lines, 444077 and 4215 are weak 
and blended. Neither of the Ba ii resonance lines, 444554 and 
4934 show more than a possibility of absorption near their posi¬ 
tions. This need not imply an underabundance of Ba or Sr as we 
anticipate significant double ionization of these two elements. 
Zr II may be weakly present, but there is no evidence for Y ii. 

6. A model using Cloudy 
6.1. Details 

It is clear that the conditions in the shell cannot be well described 
by local thermodynamic equilibrium (LTE). We have therefore 
used the versatile photoionization code Cloudy dFerland et al.l 
l2Q13ll to broaden our analysis. Ideally, we would run cloudy 
as a subroutine, and try to minimize deviations of its predic¬ 
tions against our measurements. Such an approach is beyond the 
scope of the present work. We have limited ourselves to trial and 
error comparisons of selected measureme nts or derived quan¬ 
tities, assuming d ef ault solar abundance s (lAsolund et al.l[2()()^ 
IScott et al.ll2ni5alli^ ln revesse et 1111201.511. 

The stellar background is from a iKuruczl (Il991l) Atlas9 
model. Details, and some relevant quantities are collected in Ta- 
ble|2] Properties for the = 14 000 K main sequenc e star were 
chosen with the help of data from iTorres et al.l (1201 ll) . We report 
results here for a cylindrical geometry described in the docu¬ 
mentation for Cloudy (Hazy-1, Sect. 9.6). Similar results were 
obtained with a spherical geometry with the same parameters 
apart from the cylinder height. 

It is possible to obtain column densities for any one of 371 in¬ 
dividual levels of Fe^. Our models used 256 levels, to save com¬ 
puting time. We may obtain an “excitation temperature” from 
these column densities that should be comparable to the values 
derived by the empirical curves of growth from Ti ii, Fe i, and 
Fe II. This is shown in Fig. fT3] 
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Fig. 13. Theoretical “excitation temperatures” from column densities 
for Fe II levels of the Cloudy models. The quantity g is the statistical 
weight, 27 + 1, of the levels. The larger symbols are for levels that 
correspond to the lines used in Fig. [12] The heavy line is a least squares 
solution giving T = 11134 K. Smaller symbols are for all 256 levels used 
in the model; the lighter line corresponds to T = 7026 K. 


Recall that the empirical results for Feii and Tin were 
~9164K, while that for Fei was 7750K. Theoretical and empir¬ 
ical results are summarized in Table 0 The “empirical” entries 
for neutral and ionized hydrogen are described as semi-empirical 
and shown inside brackets. While the IT formula allows us to es¬ 
timate a mean electron density (cm^^), we can only infer a col¬ 
umn density indirectly. The entry for in TablejSjcomes from 
taking the product of the IT electron density and multiplying by 
the physical thickness of the model, to obtain a column density. 
We assume the column density of protons is reasonably approx¬ 
imated by the electron column density thus obtained. 

We have an empirical estimate of the column density only 
for the n = 3 level of neutral hydrogen. We need a temperature 
to convert this to a column density for all levels of hydrogen, us¬ 
ing the Boltzmann formula. There is no obvious empirical value 
to use; indeed, because of the non-LTE conditions the meaning 
of any temperature so chosen is strictly limited to the Boltzmann 
formula and the relevant level populations. If we simply define 
a Tex as that value giving the column densities A^ 3 (H) and A^(H), 
the adopted CLOUDY model gives SOTO/Tex = 0.69. Recall the 
empirical values of 9 from Fe i and Fe ii were 0.65 and 0.55. The 
entry for H** in Table |8] was obtained using 9 — 0.69, and the em¬ 
pirical log(A? 3 (H)) = 15.3. 

Let us assume (Hanuschik 1996) that Be and shell stars have 
“disks that are intrinsically identical.” Then our basic parameters 
from Table m should be comparable with those estimated for Be 
stars. These stars display a variety of characteristics, summarized 
by Slettebak and Smith (2000). The stellar types range from late 
O (X Per) to late B iJ3 CMi). The disks are typically 5 to 20 stellar 
radii, with temperatures of 10 OOOK, and electron densities of 
10'° to lO'^ per cm. Our values for T, and radii agree well 
with these values. 

6.2. Comments 

The dominant species, e.g. Le^ are predicted reasonably well by 
the model. Lor column and particle densities in Table |8| apart 
from neutral iron, the rms deviations of the empirical from the 
9 model values is 0.33 dex, or a factor of 2.1. Neutral iron is 
an order of magnitude or more stronger than the predictions. If 
we include neutral iron, the rms deviation is 0.42 dex. A cooler 
central star could reduce this discrepancy, but the broad He i ab¬ 
sorptions would then be difficult to explain. It does not help the 


Table 8. Comparison of various empirical results with corresponding 
quantities from the model using solar abundances. Numerical values 
without units are logarithms of column densities (cm“^). Neutral species 
are indicated by a superscript ‘O’. Parenthesized numbers indicate the 
number of lines used. Values in square brackets are semi-empirical (see 
text). 


Species 

Model 

Empirical 

H« 

23.30 

[22.8] 

H+ 

23.15 

[23.5] 

N3(H) 

15.85 

15.3 

log(A'e)(cm“‘’) 

10.99 

11.2 

Sc+ 

14.54 

15.2(2) 

Ti+ 

16.43 

16.8(5) 

Cr+ 

17.11 

17.2(4) 

Fe° 

14.25 

15.4(13) 

Fe+ 

18.98 

18.8(16) 

Ni+ 

17.74 

17.6(4) 

Tex (K) 

11134 

9164 

tVAHa:)(A) 

28 

24 


ionization discrepancy to move the cloud out by a factor of 10. 
Higher or lower overall densities were also tried, with appropri¬ 
ate adjustments of the absorption path lengths. 

We note that the treatment of complex atoms of the iron 
group by Cloudy, apart from Le ii, does not involve multi-level 
models. In view of this approximate treatment, we do not con¬ 
sider the discrepancy with the Le i column density to be serious 
in the context of the present study. 

7. Conclusions 

HD 94509 is a shell star with an unusually stable shell, posing an 
extreme in the zoo of Be star spectral morphology. We have not 
been able to account for the unusual shapes of the metallic ab¬ 
sorption lines. However a simple model of the shell can predict 
the total absorption from H i, the electron density, and column 
densities of the dominant ions. Near-solar elemental abundances 
are found. The lowest Balmer lines have deep cores and symmet¬ 
rical emissions. The cores of Hy through He are virtually black, 
indicating the shell must cover the stellar disk. 
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Table A.l. Data for lines used for column densities of elements other 
than iron. References are to the oscillator strength source. 


Wavelength [A] (eV) 

log(g/) Wa (mA) 

Sen. Lawler & Dakin (1989) 

4246.82 

0.315 

0.242 

616 

4400.49 

0.605 

-0.536 

304 

4415.56 

0.595 

-0.668 

491 

4420.67 

0.618 

-l.TTi 

20 

Ti II. Pickering et al. 

(2001. 2002) 


4025.13 

0.61 

-2.14 

310 

4028.34 

1.89 

-0.96 

318 

4227.35 

1.13 

-2.24 

169 

4287.87 

1.08 

-1.82 

250 

4290.22 

1.17 

-0.85 

455 

4300.05 

1.18 

-0.44 

532 

4312.86 

1.18 

-1.10 

434 

4391.03 

1.23 

-2.28 

166 

4398.29 

1.22 

-2.78 

86 

4399.77 

1.24 

-1.19 

414 

4411.93 

1.22 

-2.52 

98 

4432.11 

1.24 

-2.81 

47 

4441.73 

1.18 

-2.27 

190 

4443.79 

1.08 

-0.72 

430 

4450.48 

1.08 

-1.52 

418 

4464.45 

1.16 

1.81 

276 

4468.51 

1.13 

-0.60 

506 

4501.27 

1.12 

-0.77 

568 

4518.33 

1.08 

-2.91 

91 

Cm. Nilsson et al. (2006) 

4195.42 

5.32 

-1.95 

41 

4242.36 

3.87 

-1.17 

324 

4558.66 

4.07 

-0.656 

616 

4588.20 

4.07 

-0.627 

545 

4634.07 

4.07 

-0.980 

454 

4812.34 

3.86 

-1.800 

341 

4824.13 

3.87 

-0.920 

500 

4836.23 

3.86 

-2.25 

137 

Niii, JCumcz (_20^ 

4015.47 

4.03 

-2.41 

58 

4067.03 

4.03 

-1.83 

269 

4244.78 

4.03 

-3.10 

17 

4362.10 

4.03 

-2.70 

59 


Appendix A: Data on iron group eiements 
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